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We present a detailed first- principles DFT study of the equation of state (EOS), energy-optimized 
geometries, phase stabilities and electronic properties of bulk crystalline CU3N, CuN and CuN2 in 
a set of twenty different structural phases. We analyzed different structural preferences for these 
three stoichiometries and determined their equilibrium structural parameters. Band-structure and 
density of states of the relatively most stable phases were carefully investigated. Further, we carried 
out GWo calculations within the random-phase approximation (RPA) to the dielectric tensor to 
investigate the optical spectra of the experimentally synthesized phase Cu3N(D0g). Obtained results 
are compared with experiment and with previous calculations. 
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I. INTRODUCTION 

In 1939, Juza and Hahn succeeded to produce CU3N 1 
for the first time 2-4 . Since then, copper nitride has been 
prepared in various techniques 3 ' 5-7 , its properties and ap- 
plications have been researched, both theoretically and 
experimentally, and it may now be considered as the most 
accessible among the noble metal nitrides 4 . 

Synthesis and reliable characterization of the prop- 
erties of a stoichiometric copper nitride constitute a 
big challenge because it is thermally unstable material 3 . 
However, this low thermal stability results in promising 
applications in optical memories and laser writing ' 6 . 

The viability of using the simple cubic stoichiometric 
CU3N films for write-once optical data storage has been 



widely explored and confirmed 8-10 , superior to other 
toxic and unstable materials in air at room temperature 
which arc used for the same purpose 9 . Also, the feasibil- 
ity of using CU3N as a coating to generate metal lines by 
maskless laser writing has been studied; where conduct- 
ing lines of a few micron in width could be generated with 
resistivities within an order of magnitude of the bulk Cu 
metal 11 . This interesting material has been suggested for 
usage in a number of nano-electronic and nano-photonic 
devices 7 . 

Depending on the total sputtering pressure and on 
the content of nitrogen gas, Hayashi et al. prepared 
four categories of sputter-deposited Cu-N films: metallic 
Cu-rich CU3N, semiconducting Cu-rich CU3N, semicon- 
ducting stoichiometric CU3N and semiconducting N-rich 
Cu3N films. In general, it has been reported that it is 
possible to achieve sub-, over- and stoichiometric copper 
nitrides 13,14 , and the effect of the nitrogen to copper ra- 
tio on the physical properties has been studied by many 
researchers 5,13-15 . 

The structural properties of CU3N in the experimen- 
tally reported cubic anti-Re03 phase are interesting on 
their own. This structure has many vacant intersti- 
tial sites like WO3. The latter could be made into a 
conductor by doping it with some metal ions 16 . This 
is very suggestive, since one may be able to engineer 
the physical properties of such technologically important 
material 15 ' 17 . In fact, the study of possible intercalated 
copper nitride alloys has been an active subject of re- 
search on its own (cf. Rcf. [5] and references therein). 

Although copper nitride possesses interesting proper- 
ties leading to different technological applications, there 
is still a large discrepancy in the formation mechanism 
and inconsistency in the experimentally reported and in 
the theoretically predicted physical properties of cop- 
per nitrides 3 ' 5 ' 7,15 ' 18 . These differences and contradic- 
tions are stemming mainly from the unstable nature (i.e. 
the metastability and low decomposition temperature) of 
copper nitride 3 ' 7,19 , the experimental conditions 15 , the 
experimental analysis methods 3 , the non-stoichiometry 
of the prepared samples 12 or the lack of knowledge of 
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the real stoichiomctry of the prepared samples 7 ; and 
from the different theoretical calculation methods and 
approximations 15 . 

Thus, the emerging potential technological applica- 
tions of copper nitride are faced by the inconsistency 
in its basic physical properties. This may explain the 
tremendously increasing interest in further studying this 
material, especially within first-principles quantum me- 
chanical approaches. Moreover, concerning its optical 
properties, only a few experiments are available in the 
literature 20 and there are very few theoretical studies. 

Motivated by all these, and searching for a wider range 
of possible applications, we present in the current work 
first-principles calculations on bulk crystalline CU3N, 
CuN and CuN 2 over a series of reported and theoreti- 
cal structural phases. The studied structural properties 
include energy- volume equation of state (EOS), equilib- 
rium lattice structural parameters, cohesive and forma- 
tion energies, relative phase stabilities, bulk modulus and 
its pressure derivative. Electronic characterization of the 
energetically most stable phases was done via the analysis 
of their band structure and their total and partial density 
of states (DOS). In order to improve the calculated elec- 
tronic structure, and to investigate the optical spectra, 
we carried out GWQ calculations within the the random- 
phase approximation (RPA) to the dielectric tensor. The 
frequency-dependent optical constants (absorption coef- 
ficient, reflectivity and refractive index spectra) of the 
experimentally reported phase Cu3N(D0g) were derived 
from the calculated frequency-dependent microscopic di- 
electric tensor. 

We hope that the present work would serve as a ref- 
erence source for meaningful comparisons which may be 
made among the largely different calculations. 



II. STOICHIOMETRIES AND CRYSTAL 
STRUCTURES 

There had been no known binary nitrides of the noble 
metals until Gregoryanz et al. 21 reported the discovery 
and characterization of crystalline PtN. From their re- 
sults, they strongly suggested that it would be possible to 
synthesize other novel nitrides with late transition metals 
such as those in the Ni and Cu groups. They claimed that 
such nitrides would have potentially intriguing physical 
properties, and that their results should stimulate further 
theoretical studies 21 . 

To the best of our knowledge, the only experimentally 
reported stoichiometries of copper nitride are CU3N 1 ' 13 ' 22 
and Cu 4 N 13 ' 22 ' 23 , while CuN and CuN 2 have not been 
observed yet. However, many transition-metal nitrides 
(TMNs) are known to form more than one nitride 24 (p. 
835). Our interest in the latter two nitride stoichiome- 
tries is based on the fact that for other late transition 
metals close to Cu in the periodic table these 1:1 and 1:2 
nitrides have been reported; as will shortly be mentioned. 

In general, the recent success in experimentally pro- 



ducing some late TMNs has motivated many researchers 
to further investigating the possibility of synthesizing 
other TMNs 25 in those reported stoichiometries and 
structures. Thus, it is of interest to know whether copper 
can form CuN and CuN 2 with the reported structures of 
other TMNs. 

In this work, we investigate CU3N in the following 
seven structures: DO3 (the fee structure of AIF93, space 
group Fm3m No. 225, Z = l) 26 , A15 (the sc structure of 
Cr 3 Si, space group Pm3n No. 223, Z = 2), D0 9 Jthc sc 
structure of the anti-Re03 (a), space group Pm3m No. 
221, Z = 1), Ll 2 (the sc structure of CU3AU, space group 
Pm3m No. 221, Z = 1), D0 2 (the bbc structure of CoAs 3 
skutterudite, space group Im3 No. 204, Z = 4), e-Fe3N 
(the hexagonal structure of e-Fe3N, space group P6322 
No. 182, Z = 2) and RhF 3 (the trigonal (rhombohedric) 
structure of FJ1F3, space group R3c No. 167, Z = 2). 

CuN was investigated in the following nine structures: 
Bl (the fee structure of NaCl, space group Fm3m No. 
225, Z = 1) which is the most popular structure for 
TMNs 27,28 and many early TMNs have been synthesized 
in this structure 29 , B2 (the sc structure of CsCl, space 
group Pm3m No. 221, Z = 1), B3 (the fee structure 
of ZnS zincblendc, space group F43m No. 216, Z = 1) 
which is the structure of the first synthesized binary ni- 
tride of the noble metals group PtN 21 , B81 (the hexag- 
onal structure of NiAs, space group P63/mmc No. 194, 
Z = 2), Bk (the hexagonal structure of BN, space group 
P63/mmc No. 194, Z = 2), Bh (the hexagonal structure 
of WC, space group P6m2 No. 187, Z = 1), B4 (the 
hexagonal structure of ZnS wurtzite, space group P63mc 
No. 186, Z = 2), B17 (the s tetragonal structure of PtS 
coopcritc, space group P4 2 /mmc No. 131, Z = 2) which 
was theoretically predicted to be the ground-state struc- 
ture of PtN 30 , and B24 (the fc orthorhombic structure of 
TIF, space group Fmmm No. 69, Z = 1). 

For CuN 2 , the following four structures were consid- 
ered: CI (the fee structure of CaF 2 fluorite, space group 
Fm3m No. 225, Z = 1) in which PtN 2 was theoreti- 
cally predicted to be stabilized 31 , C2 (the sc structure 
of FeS 2 pyrite, space group Pa3 No. 205, Z = 4), C18 
(the s orthorhombic structure of FeS 2 marcasite, space 
group Pimm No. 58, Z = 2) and CoSb 2 (the s monocline 
structure of CoSb 2 , space group P2i/c No. 14, Z = 4). 

Apart from the synthesized Cu3N(D0g), we consider 
these stoichiometries and most of these structures be- 
cause there have been some experimental reports or pre- 
vious theoretical investigations on copper nitrides or on 
other TMNs, as will be pointed out in the text. Neverthe- 
less, our first aim is to examine some possible structures 
of the infinite parameter space in order to identify the 
most energetically favorable candidates. 



III. ELECTRONIC OPTIMIZATION DETAILS 

Our electronic structure calculations were based on 
spin density functional theory (SDFT) 32,33 as implc- 
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merited in the all-electron Vienna ab initio Simulation 
Package (VASP) 34 ~ 39 . SDFT is the most widely used 40 
form of density functional theory (DFT) 41 ' 42 . Its main 
advantage over the original Hohenberg-Kohn-Sham DFT 
formalism is that it enables us to build in more of the 
actual physics into the approximate functionals with 
greater flexibility 43 . 

In solving the self-consistent Kohn-Sham (KS) 
Schrodinger-like eigenvalue equations 42 



2m„ 



dr'^\ + V ext {r) 



r - r 



(1) 



+Vj*[n(p)]Uf' k (r)= e f' k <' k (r), 



VASP expands the pseudo part of the Kohn-Sham one- 
electron spin orbitals ^' k (r) on a basis set of plane- waves 
(PWs). We included only those PWs with kinetic en- 
ergy T^Hk + G| < E cut , such that the change in total 
electronic energy and in the so-called Fermi energy Ep 
corresponding to an increase in this energy cut-off E cut 
by 100 eV is less than 3 me V/ atom and 2 meV, respec- 
tively. This is always met by E cut = 600 eV for the 
systems investigated. 

The lattice translation symmetry of a crystalline solid 
manifests itself in the quantum number k. In principle, 
there is one equation of type (1) for each k-point in the 
first Brillouin zone (BZ), and the expectation value (O) 
of any operator O is obtained by integrating its matrix 

elements k (r)|0|Vf' k (r) ) 

over all occupied bands in 

the k-space. For example, the electronic spin density 
n CT (r), which also couples equations (1) above, is given 
by 



keBZ 



N K „ 

E 



/r k ic k (r)i 5 



(2) 



where the index a indicates the spin component (j" or 
I), ^Vk.cr is the number of the occupied single-electron 
cigenstatcs ip?' (r) (with spin projection a) at each k- 
point of the sampled BZ, and f?' are the corresponding 
occupation numbers. The weights Wk should satisfy 
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(3) 



where N a is the total number of electrons with spin a. 

For performing BZ integrations, our BZs were sampled 
using r-centcred Monkhorst-Pack meshes 44 . We found 
that a 17 x 17x17 mesh corresponds to a number of k- 
points in the irreducible wedge of the Brillouin zone (IBZ) 
which is always dense enough such that any increase in 
the density of the mesh produces a change in the total 
energy less than 2 meV/atom, accompanied by a change 
in E F of less than 0.02 eV. 

For static calculations of the total electronic energy 
and the density of states (DOS), partial occupancies 



f?' were set using the tetrahedron method with Blochl 
corrections 45-47 ; while in the ionic relaxation, the smear- 
ing method of Methfessel-Paxton (MP) 48 was used. In 
doing this, the Fermi surface has been carefully treated 
and the smearing width was chosen such that the ficti- 
tious entropy - introduced by the smearing occupation 
scheme - has been kept always below 1 meV / atom. 

The generalized gradient approximation (GGA) 49 ~ 51 
was used for the exchange-correlation potentials 
V^p[?i(r)], where the Perdew-Burke-Ernzerhof (PBE) 
parametrization 52-54 is applied. The electron-ion inter- 
actions - the third term in Eq.(l) - were described by 
the projector augmented wave (PAW) method 39,55 . The 
PAW potential explicitly treats the 5 electrons of 2s 2 2p 3 
as valence electrons for nitrogen; while the 11 electrons 
of 3d 10 Ap 1 are treated as valence electrons in the Cu case, 
assuming completely filled d shell and placing all other 
electrons into the atomic core. Concerning relativistic 
effects, VASP performs a fully relativistic calculation for 
the core electrons, while for valence electrons only scalar 
kinematic relativistic effects are incorporated in the PAW 
potential via mass- velocity and Darwin corrections 38 . We 
have not considered spin-orbit interaction (SOI) of the 
valence electrons. 

The implemented blocked Davidson iteration scheme 56 
was chosen for the relaxation of the electronic degrees 
of freedom. Convergence was considered to be achieved 
when the change in the total energy and in the eigenval- 
ues between two successive self-consistent (SC) steps are 
both smaller than 1 x 10~ 4 eV. 



IV. GEOMETRY RELAXATION AND 
RELATIVE STABILITIES 

To study the energy-volume E(V) equation of state 
(EOS), and to determine the equilbrium parameters 
of each structure, we make isotropic variation of the 
cell volume while ions with free internal parameters 
are allowed to search for local minima on the Born- 
Oppenhcimer potential hyper-surface 57 , following the im- 
plemented conjugate-gradient (CG) algorithm 58 , untill 
all Hellmann-Feynman force components 59 on each ion 
are smaller than 1 x 10~ 2 eV/A. 

Cohesive energy E co h of a solid is defined relative to 
a state with all atoms neutral and infinitely separated 60 . 
Thus, in practice, it corresponds to the difference be- 
tween the crystal energy per unit cell and the total energy 
of the isolated atoms 5 



E C oh = -©(crystal) — ©(atoms). 



(4) 



Thus, cohesive energy per atom can be expressed, in our 

„61 



case, as 



E 



Cu m N„ 
coh 



(5) 



Z x (to + n) 

where Z is the number of Cu m N„ per unit cell, -E^" m 
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and E^ tom are the atomic energies, and m, n = 1, 2 or 3 
are the stoichiometric weights. 

Both crystal and atomic energies must be calculated 
at the same level of accuracy 5 ' 62 . VASP, however, cal- 
culates cohesive energies with respect to spherical non 
spin-polarised reference atoms 63 . Moreover, Vasp, in 
principle, allows only for the use of periodic systems. 
Thus, after being placed in an orthorhombic cell with 
13Axl4Axl5A dimensions, the energy of each iso- 
lated spin polarised pseudo-atom (with the same forego- 
ing electronic configuration) was calculated. The large 
dimensions of the cell ensures that there is no significant 
interaction between the atom and its images; while the 
physically incorrect spherical states are avoided by means 
of the orthorhombic symmetry (cf. Ref. [64] and Ref. 28 
therein.) 65 . T point and Gaussian smearing method with 
a small width of 0.002 eV were used, and the obtained 
atomic energies were subtracted manually from the co- 
hesive energies 

S Sud N " calculated by VASP. 
E co h represents the energy needed to decompose the 
solid into its atomic constituents 66 . Hence, those phases 
with the lower E co h are the relatively most stable. So, in 
order to investigate the relative stabilities of the phases 
under consideration, the obtained E co h as a function of 
volume V per atom were fitted 67 to a Birch-Murnaghan 
3rd-order equation of state (EOS) 68 



E{V) = E 



9V B 
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^ -1 



6-4 ^ 



^ -1 
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where Vq, Eq, Bq and B' are the equilibrium volume, 
the equilibrium cohesive energy, the equilibrium bulk 
modulus and its pressure derivative, respectively. These 
four equilibrium fitting parameters were determined by 
a least-squares method. 

Cohesive energy versus atomic volume data for the dif- 
ferent phases of CU3N, CuN2 and CuN are visualized 
graphically in Fig. 1, Fig. 2 and Fig. 3, respectively. The 
corresponding obtained equilibrium structural parame- 
ters and energetic and elastic properties are presented in 
Table I. In this table, phases are first grouped according 
to the nitrogen content, starting with the stoichiometry 
with the lowest nitrogen content CU3N, followed by the 
1:1 phases and ending with the nitrogen- richest CuN2 
ones. Within each group, phases are ordered accord- 
ing to their structural symmetry, starting from the high- 
est symmetry (i.e. space group) to the least symmetry. 
Our results are compared with available experiment and 
with previous theoretical calculations; with the calcula- 
tion methods and XC functional pointed out in the Table 
footnotes whenever appropriate. 

To deeper analyze and to compare the obtained equi- 
librium properties of the three stoichiomctries series 
with respect to one another, these quantities are de- 
picted/visualized again in Fig. 4. All quantities in this 
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FIG. 1. (Color online.) Cohesive energy _E co h(eV7atom) ver- 
sus atomic volume V (A 3 /atom) for CU3N in seven different 
structural phases. 
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FIG. 2. (Color online.) Cohesive energy _E co h(eV/atom) ver- 
sus atomic volume V (A 3 /atom) for CuN in nine different 
structural phases. 



figure are given relative to the corresponding ones of the 
fee crystalline elemental copper given in Table I. This 
will allow us to study the effect of nitridation on pure 
crystalline Cu 87 . 



A. Relative Stability: Cohesive Energy 

Considering E co h in the CU3N series, one can use Fig. 
4 to group these phases into two groups: a lower en- 
ergy (more stable than the elemental Cu) group, contain- 
ing DO9, DO2, RI1F3 and FesN structures; and a higher 
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TABLE I. Calculated and experimental zero-pressure properties of the twenty studied phases of CU3N, CuN and CUN2: Lattice 
constants (a(A), 6(A), c(A), a(°) and /3(°)), equilibrium atomic volume Vb(A 3 /atom), cohesive energy _E co h(eV^/atom), bulk 
modulus Bo(GPa) and its pressure derivative B' , and formation energy E f (eV / atom) . The presented data are of the current 
work (Pres.), experimentally reported (Exp.) and of previous calculations (Comp.). 
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a Rcf. [69]: Information is given at RTP. 
b Rcf. [70]. 

c Rcf. [71]: This is an average of 66 experimental values, at 20° C. 

d Ref. [72]: using LAP W-TB. 

e Ref. [72]: using LAPW-LDA. 

f Ref. [73]. 

s See Ref. 15 in [74]. 

h Ref. [74]: using APW-MT-LDA. 

1 Ref. [75]: Cohesive energies are given at A" and 1 atm = 0.00010 GPa; while bulk mudulii are given at room temperature, 
j Ref. [76]: using LDA 
k Ref. [76]: using BP-GGA. 
1 Ref. [76]: using PW-GGA. 
m Ref. [77]: using PAW-LDA. 
n Ref. [77]: using PAW-PW91. 
Ref. [77]: using PAW-GGA(PBE). 
p Ref. (25) in [78]: at room temperature. 
°- See Refs. (8)-(ll) in [78]. 

r Ref. [78]: using the so-called method of transition metal pseudopotcntial theory; a modified form of a method proposed by Wills and 

Harrison to represent the effective interatomic interaction. 
s Ref. [78]: using a semiempirical estimate based on the calculation of the slope of the shock velocity vs. particle velocity curves 

obtained from the dynamic high-pressure experiments. The given values are estimated at ~ 298 K . 
* Ref. [78]: using a semiempirical method in which the experimental static P — V data are fitted to an EOS form where Bg and B' are 

adjustable parameters. The given values are estimated at ~ 298 K. 
u Values obtained in the experimental work by Gallardo-Vega and Cruz 13 are between 3.810 Aand 3.830 A. 
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Ref. 
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Ref. 


[25]: 


using FP-LAPW+lo method within GGA(PBE). 
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Ref. 
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using FLAPW-LDA. 
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Ref. 
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using UPP-GGA(PBE). 
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Ref. 


[25]: 


using FP-LAPW+lo method within LDA. 
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Ref. 
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using full-potential linear muffin-tin orbital (FP-LMTO) method within GGA(PBE) 
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FIG. 3. (Color online.) Cohesive energy E co h(eV / atom) ver- 
sus atomic volume V (A 3 /atom) for CuN2 in four different 
structural phases. 



energy (significantly less stable than the elemental Cu) 
group containing Ll 2 , DO3 and A15 structures. The dif- 
ference in cohesive energy between the least stable phase 
in the lower group (FeaN) and the relatively most stable 
phase in the higher group (LI2) is 0.557 eV/atom, as one 
can see from Table 4. It is interesting to point out here 
that, except Cu3N(Fc3N), all phases in the first group are 
insulators, while all the less stable phases are metallic 88 . 

Although this simple cubic DOg phase was indeed 
found to be stable, yet one interesting result we ob- 
tained is that, in their equilibrium, the Cu3N(RhF3) 
and the Cu3N(D02) phases would be 0.001 eV/atom and 
0.002 eV/atom more stable than Cu3N(D0g), respec- 
tively. Moreover, Fig. 1 shows clearly that the E(V) 
relations of CU3N in DOg, DO2 and in RhF3 structures 
are almost identical. This marginal 89 difference in energy 
(Table I) and the almost identical E(V) curves (Fig. 1) 
may indicate the possiblity of the co-existence of these 
phases during the copper nitride synthesis process. How- 
ever, this behaviour in the EOS could be attributed to 
the structural relationships between these three struc- 
tures that being discussed in Ref. [24, p. 265] and in 
Ref. [90]. 

The bcc skutterudite structure (DO2) can be derived 
from the more symmetric sc DOg structure by simply dis- 
placing four of the N atoms located on parallel edges of 
the Cu cube to its center. This is done for two adjacent 
Cu cubes but in two vertical displacement directions, as 
nicely explained in Ref. [90]. On the other hand, to see 



the relation between DO9 and RI1F3, it is better to think 
of DOg as built of CugN octahedra (cf. Fig. I in Ref. 
[81]). A simple rotation of 60° of an octaheron about a 
shared vertex with another octahedron brings the system 
to a structure in which Cu atoms are in hep positions. 
Interested readers are referred to Ref. [24, p. 266] for 
more details. Thus, both DO2 and RI1F3 can simply be 
derived from DOg. Hence it is not surprising that these 
structural relations reflect in their EOS's and in other 
physical properties. 

The odd behaviour of the EOS of Fc 3 N with the 
existence of two minima (Fig. 1) shows that the 
first minima (to the left) is a metastable local min- 
imum that cannot be maintained as the system is 
decompressed. Cu ions are in the 6g Wyckoff positions: 
(x, 0, 0), (0, x, 0), (-x, -x, 0), (-x, 0, i), (0, -x, |), (x, x, |); 
with x = j to the left of the potential barrier (rep- 
resented by the sharp peak in Fig. 1), and x = 5 to 
the right of the peak. It may be relevant to mention 
here that Wang and Xue 86 obtained an additional local 
minimum at high pressure (lower volume) in the E(V) 
EOS of CuN(Bl). 

In the CuN series, all phases show less binding than 
the Cu(fcc), except that the simple tetragonal struc- 
ture of cooperitc (B17) is slightly more stable, with 
0.043 eV/atom lower E co h- This structure, B17, was the- 
oretically predicted to be the ground-state structure of 
PtN 30 . 

In the CuN2 nitrogen-richest phase series, we can see 
from Table I and from Fig. 4 that the phases of this 
group arc significantly more stable than all the studied 
phases, except CI, which, in contrast, is the least stable 
among the twenty studied phases. 

Comparing the relative stability of CU3N, CuN and 
CuN2, we find from Table I and from Fig. 4 that 
CuN2(C18) is the most energetically stable phase with 
0.526 eVyatom lower than the experimentally reported 
Cu 3 N(D0 9 ) phase. 



B. Volume per Atom and Lattice Parameters 

The equilibrium volume per atom Vq is an average 
quantity referring to the volume of the unit cell divided 
by the number of all atoms in the unit cell regardless of 
the type of the contained atoms. Vq is a quantity that 
is being used frequently in the literature in the calcula- 
tions from the EOS and to compare the packing of the 
different considered phases, since Vb is the inverse of the 
so-called number density. Our obtained numerical values 
are given in Table I and visualized in Fig. 4. Relative to 
the Cu(fcc) , all phases tend to slightly lower the Vb values 
except CuN(Bk) and the semiconducting CU3N phases. 

To study the structural effect of the nitrogen on the 
hosting Cu lattice, we, instead of using the commonly 
used average Vb, introduce the volume per metal atom 
Vq. In the case of CuN and CuN 2 it is numerically 
equivalent to the volume per formula unit, while for CU3N 
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FIG. 4. (Color online.) Calculated equilibrium properties of the twenty studied phases of copper nitrides. All quantities are 
given relative to the corresponding ones of the fee crystalline elemental copper given in the first row of Table I. 



8 



it equals to (volume per formula unit)/3. Hence, this 
quantity (V Cu ) may be considered as a direct measure 
of the Cu-Cu bond length and, thus, as an indicator of 
the effect of nitridation on the mechanical properties of 
the elemental Cu. That is, for a given cohesive energy, 
an increase in Vq u may/should lead to a decrease in Bq 
and vise versa, as will be seen when we discuss the trends 
in Bq values. 

In the same sub-window as Vq, obtained V Cu values 
are depicted relative to the Cu(fcc) in Fig. 4. Having a 
look at this figure, one can see a general behaviour: V Cu 
tends to increase with the increase in the nitrogen content 
and with the decrease in the structural symmetry. There 
is only one phase which has a clear odd bevaviour, that 
is CuN(Bk). It is worth to mention here that this Bk 
is not an hep structure, and we have not optimized its 
c/a ratio. Thus, this is the most open phase among all 
the investigated set. Nevertheless, all phases show an 
increase in V relative to the elemental Cu, and thus 
Cu-Cu bond is longer in all these nitrides than in the 
elemental Cu. This cannot be seen directly from the Vq 
values given in Tabcl I. 



Bulk Modulus and its Pressure Derivative 



Beside Eq and Vq, the equilibrium bulk modulus 
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appear in Eq. 6 as fitting parameters. The last parts to 
the right in Eq. 7 and Eq. 8 show that the only DFT 
calculated quantity is the total energy E (or, equivalcntly 
E co h), and that Bq and Bq are a second- and a third- 
order energy derivative, respectively. Thus, Bq and Bq 
are directly related to the curvature of the shown E(V) 
curves. 

Eq. 7 and Eq. 8 also tell us that if all phases have 
the same E co h, the Bq curve in Fig. 4 would become a 
mirror reflection-like with respect to the V Cu curve, and 
vise versa. In fact, such a general trend can be seen in 
Fig. 4. 

Compared to the parent Cu(fcc), the CuN phases tend 
to increase Bq. Such a conclusion has also been arrived 
at by Shimizu, Shirai and Suzuki 29 who calculated Bq 
for a series of 1:1 TMNs, including CuN. On the other 
hand, the considered CuN 2 phases are all, except CI, 
more compressible than the Cu(fcc). Considering the 1:3 
phases, one can easily see that the trend in E co h manifests 



itself again and divides this series into two groups: a 
group of more compressible semiconductors containing 
DOg, DO2, RhF3 and Fe3N; and a group with almost no 
change in the Cu(fcc) bulk mudulus containing LI2, DO3 
and A15 metallic phases. 

Having a look at Fig. 4, one may argue that, rela- 
tive to Cu(fcc), the lower Vq and the lower E co h of the 
CuN 2 (C2), CuN 2 (C18) and CuN 2 (CoSb 2 ) phases must 
have led to higher Bq values. Since this is not the case, 
we turn to our introduced V Cu : Fig. 4 tells us that all 
these three nitrogen-rich phases have higher V Cu rela- 
tive to Cu(fcc). Hence, V Q Cu won the competition with 
their relatively lower E co h, leading to lower Bo- This, 
again, makes sense and justifies our introduction of V Cu 
when dealing with such nitrides. 

Take, for example CuN 2 (C18). One can notice from 
Table I that, relative to Cu 3 N(D0 9 ), CuN 2 (C18) has 
about 30% less atomic volume Vq and about 14% less 
E co h, but resulting in about 30% less bulk modulus 
than Cu3N(D0g). However, if we consider our intro- 
duced V Cu , rather than the commonly used Vq, one 
can see that CuN 2 (C18) has about 57% more V Cu than 
Cu3N(D0g), which explains the lesser Bq value. In 
fact, CuN 2 (C2), CuN 2 (C18) and CuN 2 (CoSb 2 ) share al- 
most the same features when compared to Cu3N(D0g), 
Cu 3 N(D0 2 ) and Cu 3 N(RhF 3 ) as CuN 2 (C18) when com- 
pared to Cu 3 N(D0 9 ). 

Compared to the other CuN 2 phases, the relatively 
greater value of E co h of CI is overtaken by the relatively 
less V Cu value, resulting in a greater Bq than all the three 
other CuN 2 phases. It is also worth to notice from Table I 
that all the CuN phases, except the open Bk, have higher 
bulk moduli than all CU3N and CuN 2 phases, except the 
least compressible phase, CI. 

Hence, the isotropic elastic properties depend on, and 
are more appropriately described by, V Q Cu rather than Vq. 
Fig. 4 tells that the more the nitrogen content, the longer 
the Cu-Cu bond length. Thus nitrogen tends to open 
the hosting Cu lattice and to reduce the bulk mudulus. 
Nevertheless, bulk mudulus is a result of the competition 
between the change in V Q Cu and the change in E co h- 

Physically, the pressure dependence of Bq can be quan- 
tified via its pressure derivative B' given by Eq. 8. Ex- 
cept for the last three CuN 2 phases in Fig. 4, all phases 
show almost equal sensitivity. However, the change 
in Bq of the elemental Cu would be greater than all 
these phases. On the other hand CuN 2 (C2), CuN 2 (C18) 
and CuN 2 (CoSb 2 ) show high elastic sensitivity to any 
isotropic pressure. It is worth to mention here that this 
quantity is a measurable quantity 78 , but we cannot find 
any experimental value for the synthesized Cu3N(D0g) 
phase. 

Table I shows that our obtained Bq and Bq values agree 
well with many other theoretical works. However, clear 
differences between the current and, and among, the the- 
oretically obtained values exist. This can be traced back 
to two factors: 

• From Eq. 7 and Eq. 8, one needs only to calculate 
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E co h from DFT. Hence, B' values depends on the 
DFT algorithm/method and functional. For exam- 
ple, in Table I, the GGA (e.g. Refs. 25 and 85) 
calculated Bo values of Bl, B2 and B3 are all lower 
than the obtained LDA (e.g. Refs. 25 and 29) 
corresponding ones. This is, in fact, a well-known 
feature of GGA compared to LDA. That is, relative 
to the latter, the former undcrstimatcs the cohesion 
and thus ovcrstimates Vq resulting in a lower Bq. 

• Yet, there are considerable differences among the 
GGA values and among the LDA values! Recalling 
that Bo and B' are obtained from the EOS fitting, 
the accuracy in the obtained values depends not 
only on the accuracy of the DFT calculated E co h, 
but also on the numerical fitting; e.g. number of 
points around Vq and how far these points are from 
both Vb and E$. Moreover, B' Q is numerically more 
sensitive than Bq 64 , that why, we believe, there is 
no even clear/general trend/systematic in the cal- 
culated B„ values of Table I. 



D. Relative Stability: Formation Energy 

Since the driving force for the formation of a solid 
may lie in the energy released when the elements con- 
dense to form the solid 91 , it has become a common prac- 
tice in ab initio studies to calculate the so-called en- 
ergy of formation in order to test the formation pos- 
sibility of materials under consideration 92 as well as a 
measure/indicator of the relative stabilities of the phases 
under consideration 5,93 . 

Within first-principles calculations, the formation en- 
ergy Ef can be calculated from the difference of the cohe- 
sive energies (or enthalpies) of products i? co h (products) 
and reactants i? co h (reactants) 94 ' 95 

Ef = 2^ -Ecoh (products) — -Ecoh (reactants); (9) 

where -E C oh (reactants) should be given in the reference 
states of the reactants; that is, at their most stable state 
at specified temperature and pressure 94 . 

In our case, if we assume that our product Cu m N n 
results from the interaction between the gaseous diatomic 
molecular N2 and the solid Cu metal in its reference fee 
Al structure through/via the chemical reaction 

TO Cu solid + ^Nf s — > Cu rn C lid , (10) 

then Eq. 9 above can be written for formation energy 
per atom as (cf. Refs. 5, 93, and 96): 

£/(Cu m Nf d ) = £ coh (Cu m Nf id ) 

m£ coh (Cu solid ) + f £ coh (Nf as ) (11) 
to + n 



where m,n ~ 1,2,3 are the stoichiometric weights and 
-E'coh(Cu m N^ olld ) is, again, the cohesive energy per atom 
calculated via Eq. 5 97 . 

To determine the equilibrium cohesive energy of the 
elemental metallic copper £coh(Cu solid ) in its well-known 
fee Al structure (space group Fm3m No. 225) 73,98,99 , we 
followed the same procedures described in Sec. Ill and 
Sec. IV. The obtained structural and cohesive properties 
of the bulk Cu are placed in the first row of Table I, where 
they show excellent agreement with experiment and good 
agreement with the theoretically predicted ones. 

The total energy of the gaseous diatomic molecular ni- 
trogen (i? co h(N| as )), was obtained by placing one N atom 
at a corner of a cubic cell with a = 14 A, while the second 
atom is displaced along the diagonal direction and was 
allowed to move during the relaxation. E cut of 800 eV, 
r point and Gaussian smearing method with width of 
0.05 eV were used. Molecular and cohesive energies were 
calculated in the same way as described in Sec. IV. Ob- 
tained cohesive and structural results are given in Table 
II and compared with experiment and with previous cal- 
culations therein. 

TABLE II. Calculated and experimental cohesive energy 
E co h_{eV) and bond length d(A) of the gaseous diatomic 
molecular nitrogen (Nf as ). The presented data are of the 
current work (Pres.), experimentally reported (Expt.) and of 
previous calculations (Comp.). 

Pres. Expt. Comp. 

E cob {eV) -10.392 -(9.797658 ± 0.0061149) a 10.49°, 11.75 c , 11.57", 10.69 e 
d(A) 1.113 (1.0976 ± 0.0002) f 1.102 b , 1.085 c , 1.095 d , 1.095° 

a This bond strength in nitrogen diatomic molecule is taken from 
Ref. [70, p. 9:55] where it is given there as 
(945.33 ± 0.59 KJ/mol) with positive sign convention and at 
298 K. In p. 9:76 of the same reference, the force constant for 
bond stretching in nitrogen diatomic molecule is given to be 
22.95 N/cm. 

b Ref. [100], PBE(GGA)-LAPW, with spherical ground-state 

density of the N free atom. 
c Ref. [100], LDA-PP, with spherical ground-state density of the 

N free atom. 

d Ref. [100], LDA-LAPW, with spherical ground-state density of 

the N free atom. 
° Ref. [100], PBE(GGA)-PP, with spherical ground-state density 

of the N free atom. 
f Ref. [101, p. S7]. 

The obtained formation energies Ef of the twenty re- 
laxed phases are given in Table I and shown graphically in 
Fig. 4. All these values are positive; which means that all 
these twenty phases are, in principle, thermodynamically 
unstable (endothermic) . However, these results have to 
be interpreted with some caution: 

• Many other theoretical calculations found posi- 
tive formation energy for experimentally synthe- 
sized transition metal nitrides; e.g. OSN2 102 (us- 
ing PP and PBE-GGA), PtN 2 103 (using PP and 
PBE-GGA) and InN 100 (using PP with LDA and 
different GGAs) 104 (using PW91-GGA). 
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• The fact that we obtained a positive formation en- 
ergy for the successfully synthesized Cu 3 N(D0 9 ) 
phase means that it may be possible that other 
phases can still be synthesized, and it may indi- 
cate that there is a problem with our calculations 
method (i.e. with the approximations) and/or with 
the physical conditions assumed for the calculations 
(i.e. pressure and temperature); see below. 

• These positive values are the result of static DFT 
calculations (T = K) at equilbrium volume 
(P = GPa); while the fact is that most, if not 
all, of the successfully synthesized TMNs were ob- 
tained by subjecting their parent elements to ex- 
treme conditions of pressure and temperature (c.f. 
Rcf. [102]). 

• Referring to Table II, the difference between our 
calculated cohesive energy of N 2 and experiment 
(Ref. [70]) is about -0.297 eU/atom; while for 
the bulk Cu the difference is about 0.016 eU/atom. 
Now, using Eq. 11 with this significant overes- 
timation of S co h(N| as ) and the reasonable value 
of -E C oh(Cu solld ) will result in underestimation of 
.E/(Cu m N^ olld ). This contribution has to be consid- 
ered, as an apparent shortcoming of the PBE-GGA, 
whenever one deals with a dimeric crystal 64,100 . 

• Nevertheless, since all formation energies are cal- 
culated as the difference between the ab initio co- 
hesive energies, which in turn are calculated at the 
same level of accuracy, one can still use these for- 
mation energies to measure the relative thermody- 
namic stabilities of these structures. That is, the 
lower the formation energy, the lower the propen- 
sity to dissociate back into the constituent elements 
Cu and N 2 (c.f. Rcf. [102]). 

• Moreover, because Cu3N(D0g) has been synthe- 
sized, we can take it as a reference measure of sta- 
bility. It is also worth to recall here that experi- 
ment found Cu3N(D0g) to be metastable at room 
temperature 3 . 

Relative to each other, and within each series, Ef of 
the twenty phases shows almost the same trend as E co h- 
However, the CuN phases tend to be relatively less sta- 
ble than the CU3N and CuN 2 phases, except the odd CI 
phase. In fact, C2, C18 and CoSb2 are the most sta- 
ble and share almost the same features when compared 
to Cu 3 N(D0 9 ), Cu 3 N(D0 2 ) and Cu 3 N(RhF 3 ). This may 
agree well with Armenta and Soto 5 who proved, from 
the study of formation energy, that the metallic phases 
of copper nitrides would be more stable than the semi- 
conducting phase. 



E. More Comparison with Experiment and with 
Theory 

Comparing our obtained results with experiment, one 
can see from Table I that the lattice parameter a of 
Cu 3 N(D0c)) was reproduced very well. Excellent agree- 
ment with previous calculations is also clear, though, 
with respect to experiment (or: with respect to each 
other), the common overcstimation of a by GGA and 
the underestimation of a by LDA (c.f. Rcfs. 15, 79, and 
81; and Rcfs. 105 and 106) is showing up. 

Using the full-potential (linearized) augmented plane 
waves plus local orbitals (FP-LAPW+lo) method within 
LDA and within GGA, Kanoun and Said 25 studied the 
E(V) EOS for CuN in the Bl, B2, B3 and B4 struc- 
tures. While within GGA, they found equilibrium lattice 
parameters which are in excellent agreement with ours, 
their obtained LDA lattice parameter values show the 
common underestimation with respect to our and their 
GGA values (see Table I). Also, the relative stabilities 
of these phases they arrived at agree well with ours, and 
they concluded that B3 is the ground-state phase of CuN 
and is metalic. 

Shimizu, Shirai and Suzuki 29 performed first-principles 
calculations using full-potential linearized augmented- 
plane-wave (FLAPW) method in the framework of LDA 
and found that CuN(Bl) is less than 0.20 eU/atom more 
stable than CuN(B3), while we found that CuN(B3) is 
0.043 eU/atom (GGA) more stable than CuN(Bl). Some 
of their findings are shown in Table I; and, again, their 
predicted LDA lattice constants are slightly less than our 
GGA values, while their obtained bulk moduli are over- 
estimated when compared to ours. 

Using full-potential linear muffin-tin orbital (FP- 
LMTO) method within GGA(PBE), Wang et al. 86 stud- 
ied the E(V) EOS of CuN in the Bl, B2, B3, B81 and 
B4 structures. Their obtained equilibrium lattice param- 
eters and bulk muduli arc included in Table I which show 
good agreement with ours. However, Fig. 2 shows that, 
within this parameter space, equilibrium cohesive energy 
decreases as B2, B81, B4, Bl and B3. This is consistent 
with Wang et al. but B81 and B4 are swaped. Neverthe- 
less, B3 is the most stable in both works, contrary to the 
findings of Shimizu, Shirai and Suzuki 29 . 

Whatever the case, in our wider parameter space, Fig. 
2 and Table I reveal that CuN(B17) is 0.17 eV/atom 
and 0.21 eU/atom (GGA) more stable than CuN(B3) 
and CuN(Bl), respectively. It may be worth to mention 
again here that B17 was theoretically predicted to be the 
ground-state structure of PtN 30 . 

Using norm-conserving ultra-soft pscudopotcntial 
within GGA and the so-called BFGS scheme for geom- 
etry optimization, Bouayed et al. 85 studied CuN in the 
Bl and B3 structures, and CuN 2 in the CI structure. 
Although their obtained lattice constants (given in Table 
I) are in good agreement with our findings, the notica- 
blc difference in bulk moduli may be traced back to the 
numerical fitting (see Sec. IV C). 
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V. ELECTRONIC PROPERTIES 

Band structure (i.e. ef (k) curves) and spin-projected 
total and partial (i.e. orbital resolved) density of states 
(DOS) of the energetically most stable phases are pre- 
sented in Figs. 5, 6, 7, 8 and 9. Spin-projected 
total density of states (TDOS) are shown in subfig- 
ure (b) in each case. In all cases, TDOS's arc com- 
pletely symmetrical in majority and minority spins. 
That is, electrons occupy the majority and minority 
spin bands equally and result in a zero total spin mo- 
ment and a zero spin-polarization ratio: SPRdos(-E0 = 
| (D t (E) - D^E)) I (D^(E) +D i (E))\. That is why it 
was sufficient only to display spin-up partial density of 
states (PDOS) and spin-up band structures. To properly 
show details of the electronic structure of these phases, 
we plotted the energy bands along densely sampled high- 
symmetry string of neighbouring points in the k-space; 
while displaying the Cu(s,p, d) and N(s,p) partial DOS 
allows us to extract information about the orbital char- 
acter of these bands. 

Energy bands ef (k), total density of states (TDOS) 
and partial (orbital-resolved) density of states (PDOS) 
of Cu3N(D0g) are shown in Figs. 5. It is clear that 
Cu3N(D0g) presents insulating character in its spin band. 
It has its valence band maximum (VBM) at R which lies 
0.13 eV below the Fermi energy Ep, and its conduc- 
tion band minimum (CBM) at M which also lies 0.13 
eV above Ep, resulting in a narrow indirect gap of 0.26 
eV . By looking at the PDOS plots, energy bands can be 
divided into three parts: a deep band around ~ — 16 eV 
below Ep consists mainly of N(2s), a broad group of 12 
valence bands with ~ 8 eV of width that comes mostly 
from the 3d electrons of Cu plus smaller contribution 
from N(2p), and the conduction bands. 

Our obtained band structure of Cu3N(D0g) agrees 
qualitatively well with many previous theoretical 
works 2 ' 15,79 ; however, depending on the calculation 
method, the value of the indirect band gap of CusN(DOg) 
was predicted to be 0.23 eV (LAPW+LDA) 2 , 0.25 
eV (FP-LAPW+GGA(PBE)) 17 > 79 , 0.355 eV (UPP- 
GGA(Perdew-Wang)) 15 and 0.9 eV (LCAO+LDA) 2 . Al- 
though our predicted band structure and gap value agree 
well with many theoretical calculations, the experimen- 
tally reported values of the energy gap E g of Cu3N(D0g) 
are larger, as will be discussed in Sec. VI. This is a well 
known drawback of Kohn-Sham DFT-based calculations 
to understimate the band gap. 

Fig. 6 shows that the top of the valence band and 
the bottom of the conduction band of Cu3N(RhF3) are 
about to touch the Fermi level at (r, — 0.13 eV) and 
(£, 0.15 eV), respectively. Hence, we have an indirect 
band gap of 0.28 eV with Ep slightly shifted towards the 
VBM. 

The third most stable candidate in this 3:1 series is 
Cu 3 N(D02). Its band structure (shown in Fig. 7) has the 
VBM at (if, -0.14 eV), and the CBM at (r,0.14 eV), 
resulting in an indirect energy gap of 0.28 eV. 



The orbital resolved density of states (PDOS) of both 
Cu 3 N(D02) and Cu3N(RhF 3 ) share the same qualita- 
tive features with those of CusN(DOg). As expected, 
the structural relation between these three phases are 
reflected into their electronic properties. 

Energy bands ef (k), total density of states (TDOS) 
and partial (orbital-resolved) density of states (PDOS) 
of CuN(B17) are shown in Figs. 8. It is clear that 
CuN(B17) would be a true metal at its equilibrium. The 
major contribution to the very low TDOS around Fermi 
energy Ep comes from the 2p states of the N atoms. Be- 
neath Ep lies a band with ~ 8 eV of width, in which 
one can notice that the main contribution is due to the 
mixture of Cu(3g?) states with N(2p) states. The N(2s) 
states dominate the deep lowest region, while the unoc- 
cupied states stem mainly from the N(2p) states. Fermi 
surface intersects two bands: a lower one in the M-T- 
X-R directions, and a higher band in the T-M-A and 
M-X directions. Hence, Ep is not a continous surface 
contained entirely within the first BZ. 

So, in CuN(B17), the nitridation effects on the pure Cu 
can be summarized from previous sections as: significant 
increase in the volume per atom Vq, increase in the bulk 
modulus resulting in a less compressible material than 
the pure metal, while the metallic character is preserved 
. Similar results for CuN(B3) were also arrived at by 
other researchers 85 . 

It may be worth to mention here that Bl 25 ' 85 , B2 25 ' 85 , 
B3 25 and B4 25 phases of CuN were also theoretically pre- 
dicted to be metallic. 

With Ep crossing the finite TDOS, Fig. 9 shows that 
CuN 2 (C18) is metallic at K. The PDOS reveals that 
the major contribution to the TDOS at Ep comes from 
the N(2p) states with minor contribution from the Cu(3cf) 
states. Compared to CuN(B17), a new feature of this 1:2 
nitride is the contribution of N(2p) states to N(2s) states 
at the deep lowest region. However, variation in N(2s) 
energy with respect to k is smaller than the variation of 
N(2p) states, resulting in a narrower and higher PDOS. 
It may be instructive to mention here that CuN2(Cl) 
phase was also found to be metallic 85 . 

A common feature between all the studied cases is the 
higly structured, intense and narrow series of peaks in the 
TDOS valance band corresponding to superposition of 
N(2p)-states and Cu(3d)-states. In their k-space, Cu(3d) 
energies show little variation with respect to k; hence the 
Hove singularities-like sharp features. 

To summarize this section, we have found that the 
most stable phases of CuN and CuN 2 are metallic, while 
CU3N is a semiconductor. This finding agrees well with 
literature, specially with Armenta and Soto 5 who pre- 
dicted theoretically that the semiconducting state is sen- 
sitive to the nitrogen concentration and changes to metal- 
lic when the composition is out of the ideal nitrogen to 
copper ratio, x = 1/3. Armenta and Soto, who studied 
the effect of introducing N atoms in one by one basis to 
the bulk cubic Cu matrix, also pointed out that as x in- 
creases, the TDOS at Ep, due to both N and Cu atoms, 
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FIG. 5. (Color online.) DFT calculated electronic structure for CU3N in the DOg structure: (a) spin-projected total density 
of states (TDOS); (b) partial density of states (PDOS) of Cu{s,p,d) orbitals in Cu 3 N; (c) PDOS of N(s,p) orbitals in Cu 3 N, 
and (d) band structure along the high-symmetry k-points which are labeled according to Ref. [107]. Their coordinates w.r.t. 
the reciprocal lattice basis vectors are: M(0.5, 0.5, 0.0), T(0.0, 0.0, 0.0), X (0.0, 0.5, 0.0), #(0.5, 0.5, 0.5). 
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FIG. 6. (Color online.) DFT calculated electronic structure for Cu 3 N in the RI1F3 structure: (a) spin-projected total density 
of states (TDOS); (b) partial density of states (PDOS) of Cu(s,p, d) orbitals in Cu 3 N; (c) PDOS of N(s,p) orbitals in Cu 3 N, 
and (d) band structure along the high-symmetry k-points which are labeled according to Ref. [107]. Their coordinates w.r.t. 
the reciprocal lattice basis vectors are: F(0.5, 0.5, 0.0), Q(0.375, 0.625, 0.0), 5(0.5,0.75,0.25), 2(0.5,0.5,0.5), T(0.0, 0.0, 0.0), 
£(0.0,0.5,0.0), F(0.25,0.5, -.25), £(0.0, 0.5, -.5). 



increases as well. Concerning this point, our findings are 
in excellent agreement with theirs, since we found for the 
most stable phases (i.e. CU3N, CuN and CuN2, respec- 
tively) that TDOS(.t = 1/3) = 0, TDOS(a; = 1) - 0.70 
and TDOS(x = 2) ~ 0.85, in relative arbitrary units. 



VI. GWA CALCULATIONS AND OPTICAL 
PROPERTIES 

Although a qualitative agreement between DFT- 
calculated optical properties and experiment is possible, 
accurate quantitative description requires treatments be- 
yond DFT level 108 . Another approach provided by many- 
body perturbation theory (MBPT) leads to a system of 
quasi-particle (QP) equations, which can be written for 



a periodic crystal as 109 111 

Practically, the wave functions are taken from the DFT 
calculations. However, in consideration of computational 
cost, we used a less dense mesh of k-points (12 x 12x 12). 
The term E(r, r'; e^) is the self-energy which contains 
all the exchange and correlation effects, static and dy- 
namic, including those neglected in our DFT reference 
system. In the so-called GW approximation, E is given 
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FIG. 7. (Color online.) DFT calculated electronic structure for CU3N in the DO2 structure: (a) spin-projected total density 
of states (TDOS); (b) partial density of states (PDOS) of Cu(s,p,d) orbitals in Cu 3 N; (c) PDOS of N(s,p) orbitals in Cu 3 N, 
and (d) band structure along the high-symmetry k-points which are labeled according to Ref. [107]. Their coordinates w.r.t. 
the reciprocal lattice basis vectors are: T(0.0, 0.0, 0.0), iV(0.0, 0.0, 0.5), P(0. 25, 0.25, 0.25), #(0.5, -.5, 0.5). 




FIG. 8. (Color online.) DFT calculated electronic structure for CuN in the B17 structure: (a) spin-projected total density 
of states (TDOS); (b) partial density of states (PDOS) of Cu(s,p,d) orbitals in CuN; (c) PDOS of N(s,p) orbitals in CuN, 
and (d) band structure along the high-symmetry k-points which are labeled according to Ref. [107]. Their coordinates 
w.r.t. the reciprocal lattice basis vectors are: X(0.0, 0.5, 0.0), M(0.5, 0.5, 0.0), T(0.0, 0.0, 0.0), Z(0.0, 0.0, 0.5), A(0.5, 0.5, 0.5), 
#(0.0,0.5,0.5). 




FIG. 9. (Color online.) DFT calculated electronic structure for CuN2 in the C18 structure: (a) spin-projected total density 
of states (TDOS); (b) partial density of states (PDOS) of Cu(s,p, d) orbitals in CuN 2 ; (c) PDOS of N(s,p) orbitals in CuN 2 , 
and (d) band structure along the high-symmetry k-points which are labeled according to Ref. [107]. Their coordinates 
w.r.t. the reciprocal lattice basis vectors are: T(0.0, 0.0, 0.0), X(0.0, 0.5, 0.0), 5(-.5, 0.5, 0.0), Y (-.5, 0.0, 0.0), £(0.0,0.0,0.5), 
1/(0.0,0.5,0.5), #(-.5,0.5,0.5), T(-.5, 0.0, 0.5). 
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FIG. 10. (Color online.) Normal-incidence frequency- 
dependent optical constants of Cu3N(D0g) obtained using 
Eqs. 15-17 and GWo eigenvalues. 



in terms of Greens function G as 

Sgw = 3 J de'G(r, r'; e, e')W(r, r'; e), (13) 

where the screened interaction W is related to the bare 
Coulomb interaction v through 



W(r,r';e) =j J dr^ 1 (r, n; e)v(r 1} r'), 



(14) 



with e the dielectric function. We followed the GWq self- 
consistent routine on G, in which the QP eigenvalues are 
updated in the calculations of G, while W is kept at the 
DFT level. After the final iteration of G, e is calcu- 
lated, within the so-called random phase approximation 
(RPA) 112 , using the updated QP eigenvalues 63,110 ' 111 . 
From the real e re (uj) and the imaginary Ei m (uj) parts 
of this frequency-dependent microscopic dielectric ten- 
sor one can derive all the other frequency-dependent di- 
electric response functions, such as absorption coefficient 
a (uj), reflectivity R (uj) and refractive index n (uj): 



R(e>) 



, (ui) +je im (uj)}- 



1 



[e re (uj) + je im (cj)] 2 
a(u,)=V2u,([s 2 re (w)+ jeLM]* 



+ 1 



(15) 



H) 5 (16) 



'4M+4H] 5 + £re ( w ) (i7) 



Fig. 10 displays the real and the imaginary parts of 
the frequency-dependent dielectric function £rpa(cl>) of 
Cu3N(D0g) and the corresponding derived optical con- 
stants within the optical region 113 . The real part e re (uj) 
shows an upward trend before 0.3 eV, reaches a maxi- 
mum value at 1 eV and generally decreases after that. 



The imaginary part Si m (uj) has two main peaks located 
at ~ 1.3 eV and ~ 2.0 eV. Niu et al. 82 carried out 
DFT(UPP-GGA) calculations and derived an e(u) spec- 
trum with a real part that shows an upward trend before 
1 eV, reaches the maximum value at 1 eV (as ours) and 
decreases after that. The imaginary part they obtained 
has two main peaks, in the range eV to 5 eV , located at 
2.07 eV and 3.51 eV. By analyzing the DOS, Niu et al. 
claimed that these two peaks are mainly due to the elec- 
tron transition from the Cu(3e?) band to the conduction 
band. 

Refractive index n (uj) spectrum (Fig. 10) shows al- 
most the same frequency dependance as e re (uj). Compar- 
ing our obtained n (uj) with the experimental results of 
Gordillo et al. 7 , one can see a kind of general qualitative 
agreement between the experimental curve and the theo- 
retical one, represented by the increase in n with increase 
in the photon energy till reaching a maximum, then fol- 
lowed by a generally decreasing behaviour. However, the 
experimental peak of n is at ~ (1.4 eV,3.65), while our 
obtained peak is at ~ (1.2 eV, 4.5). Surprisingly, the 
DFT(UPP-GGA)-refractive index spectrum derived by 
Niu et al. 82 shows far better agreement with the experi- 
mental results of Gordillo et al. 7 and they got the peak 
at - (1.6 eV,3.9). 

From the absorption coefficient a (uj) spectrum (Fig. 
10), it can be seen that Cu3N(D0g) starts absorbing pho- 
tons with ~ 0.75 eV energy. Hence, it is clear that 
GWq calculations give a band gap of ~ 0.75 eV, which 
is a significant correction to the obtained DFT value. 
Our presented a (uj) spectrum agrees qualitatively well 
with the experimental work of Gordillo et al. 7 , who, 
in the ~ (0.6 — 1.4) eV region, obtained a smooth 
exponential-like curve. However, their obtained a (uj) 
reaches 1 x 10 5 cm" 1 in the visible range before 1.4 eV, 
while ours (Fig. 10) never reaches such a value before 

3.3 eV. Niu et al. 82 got a curve that reaches this value 
at ~ 3.6 eV. However, their a (uj) spectrum starts to be 
non-zero from ~ 0.71 eV\ 

Gordillo et al. 7 prepared nearly stoichiometric copper 
nitride polycrystalline films having nitrogen contents of 
(27 ± 2)% with lattice parameter a = 3.8621 A. They re- 
ferred to it in their article as stoichiometric CU3N, and, at 
room temperature and with orientation along the (10 0) 
crystallographic axis, they carried out some optical mea- 
surements and fitted the obtained data. From the fits, 
they managed to derive the refractive index and the ab- 
sorption coefficient, while reflectance was measured di- 
rectly. So, although our calculated optical properties 
show partial agreement with this experimental work, dis- 
crepancies may be attributed due to the lack of knowl- 
edge of the exact stoichiometry of the prepared samples. 
Wang, Nakamine, and Hayashi 12 also prepared nearly 
stoichiometric Cu 3 N films at 67 Pa. Their experimen- 
tally obtained a (uj) spectra reach 1 x 10 5 cm -1 at about 

2.4 eV. However, a leaves the zero level only at about 
(1.8 - 1.9) eV. 

On the other hand, the used E cut = 290 eV and 8x8x8 
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k-mesh in the DFT(UPP-GGA) calculations by Niu ct 
al. 82 may not be sufficient to reproduce qualitatively sim- 
ilar spectra as those we obtained from GWA calculations. 
No quantitative correspondence is to be expected. 

Experimentally reported values of the Cu3N(D0g) opti- 
cal gap spread over a wide range 3,7 . Some of these values 
are: (0.25 - 0.83) eV 14 , 1.30 eV 6 < 83 , (1.7 - 1.84) eV 80 , 
1.85 eV 4 and (1.8 - 1.9) eV 12 . Hence, although GWA 
calculations give a band gap of ~ 0.75 eV, which is a sig- 
nificant correction to the obtained DFT value of 0.26 eV, 
the difference between theory and experiment is still con- 
siderable. This is a well known problem with GWA cal- 
culations on top of DFT eigenvalues and eigenstates that 
correspond to a very small bandgap compared to exper- 
iment. 

Given that standard DFT functionals severely under- 
estimate the band gaps while the Hartree-Fock (HF) ap- 
proximation overestimates them 114 , a potential solution 
is to combine local or semilocal DFT exchange with a 
portion of nonlocal exact exchange thereby construct- 
ing the so-called hybrid functional 115 . Much improved 
band gaps can be obtained by screening the nonlocal HF- 
type portion of exchange potential with a suitable screen- 
ing parameter 116,117 . A more appropriate approach is 
to apply the partially self-consistent GW method on 
eigenvalues obtained using hybrid functionals or DFT+U 
schemes 118 which readily provide better band gaps for 
insulators 119 . 



VII. SUMMARY AND CONCLUSIONS 

DFT-based first-principles calculations on bulk crys- 
talline Cu 3 N, CuN and CuN 2 over a series of twenty 
structural phases have been successfully carried out. 
The studied structural properties include energy- volume 
equation of state (EOS), equilibrium lattice structural 
parameters, cohesive and formation energies, relative 
phase stabilities, bulk modulus and its presssure deriva- 
tive. Electronic characterization of the energetically most 
stable phases was done via the analysis of their band 
structure and their total and partial density of states 
(DOS). Further, we carried out GWo calculations within 
the random-phase approximation (RPA) to the dielec- 
tric tensor £rpa(w). The frequency-dependent optical 
constants (absorption coefficient, reflectivity and refrac- 
tive index spectra) of the experimentally reported phase 
Cu 3 N(D0g) were derived from the updated £rpa(w). Ob- 
tained results were compared with experiment and/or 
with previous calculations whenever possible. The main 
conclusions which we can derive from all these calcula- 
tions are the followings: 

• The calculated lattice constants are in good agree- 
ment with experiment and with theory. 



• From the obtained cohesive energies, the energeti- 
cally most stable phases are D0g, B17, and C18 in 
the CU3N, CuN and CuN 2 stoichiometric series, re- 
spectively. However, other CU3N phases show simi- 
lar stability to Cu3N(D0g) and may present during 
the nitridation process. 

• Including the successfully synthesized Cu3N(D0g) 
phase, all obtained formation energies are positive, 
yet they can be used to measure the relative ther- 
modynamic stabilities of these phases. 

• Although CuN2(C18) is the most bound phase, its 
tendency to decompose back into its elemental con- 
stituents is more than the less bound CU3N phases. 

• The volume dependence of the bulk modulus is 
more precisely described by the change in volume of 
the Cu sublattice rather than the common average 
atomic volume of the nitride. 

• The most stable CU3N phases are predicted to be 
indirect-gap semiconducing materials with lower 
bulk modulus than the pure metal, while CuN(B17) 
preserves the metallicity and improves the bulk 
modulus. However, the CuN 2 (C18) phase substan- 
tially increases the compressibility while preserving 
the metallicity. 

• Our GWA calculated optical properties show par- 
tial agreement with experiment and with the avail- 
able theoretical work (Ref. [82]). Discrepancies are 
probably due to the lack of knowledge of the exact 
stoichiomctry of the prepared samples, and due to 
the big difference in the used plane waves cut-off 
energy and in the density of the k-mesh. Conver- 
gence criterion with respect to these two quantities 
have not been referred to in Ref. [82]. 

• Our GWA calculated energy gap of Cu 3 N(D0g) 
shows significant improvement over the calculated 
DFT value. 

We hope that the present work would serve as a reference 
source for meaningful comparisons which may be made 
among the largely different calculations. 

ACKNOWLEDGMENTS 

All GW calculations and some DFT calculations were 
carried out using the infrastructure of the Centre for 
High Performance Computing (CHPC) in Cape Town. 
Suleiman would like to acknowledge the support he re- 
cieved from Wits, DAAD, AIMS and SUST. Many thanks 
to the ASESMA family, and special thanks to Dr Kris 
Delaney and Sinead Griffin for their invaluable help and 
useful discussions. 



16 



Corresponding author: suleiman@aims.ac.za 6, 9 (2004), ISSN 1293-2558, 

Current address: Department of Physics, University of http://www.sciencedirect.com/science/article/pii/S129325580; 

Pretoria, Pretoria, South Africa. 18 J.-N. Ding, N.-Y. Yuan, F. Li, G.-Q. Ding, 

R. Juza and H. Halm, Z.-G. Chen, X.-S. Chen, , and W. Lu, 

Zeitschrift fr anorganische und allgemeine Chemie Journal of Chemical Physics 131, 174102 (2009), ISSN 

241, 172 (1939), ISSN 1521-3749, 1089-7690, http://dx.doi.Org/10.1063/l.3257899 

http://dx.doi.org/10.1002/zaac.19392410204 19 G. Yue, P. Yan, and J. Wang, Journal of Crystal Growth 

U. Hahn and W. Weber, 274, 464 (2005), ISSN 0022-0248, 

Phys. Rev. B 53, 12684 (May 1996), http: //www . sciencedirect . com/science/article/pii/S002202480^ 

http://link.aps.org/doi/10.1103/PhysRevB.53.126 8 4 20 G. Soto, J. Daz, and W. de la Cruz, 

A. Ji, R. Huang, Y. Du, C. Li, Y. Wang, Materials Letters 57, 4130 (2003), ISSN 0167-577X, 

and Z. Cao, Journal of Crystal Growth http://www.sciencedirect.com/science/article/pii/S0167577X0; 

295, 79 (2006), ISSN 0022-0248, 21 E. Gregoryanz, C. Sanloup, M. Somayazulu, 

http://www.sciencedirect.com/science/article/pii/S0022024806BMB44 G. Fiquet, H. kwang Mao, and 

Y. Du, A. Ji, L. Ma, Y. Wang, R. J. Hemley, Nature Materials 3, 294 (2004), 

and Z. Cao, Journal of Crystal Growth http://www.nature.com/nmat/journal/v3/n5/full/nmatlll5.html 

280, 490 (2005), ISSN 00 22-0248 , 22 F. Gulo, A. Simon, J. Khler, and R. K. Kre- 

http : //www. sciencedirect . com/science/article/pii/S002202A8fi§004264 Angewandte Chemie International Edition 

M. G. Moreno-Armenta and G. Soto, 43, 2032 (2004), ISSN 1521-3773, 

Solid State Sciences 10, 573 (2008), ISSN 1293-2558, http://dx.doi.org/10.1002/anie.200353424 

http : //www. sciencedirect . com/science/article/pii/S1293 2 i538070B2Q8tier, K. Bang, and B. Giessen, 

S. Ghosh, F. Singh, D. Choudhary, Materials Science and Engineering: A 

D. Avasthi, V. Ganesan, P. Shah, and 117, LI (1989), ISSN 0921-5093, 

A. Gupta, Surface and Coatings Technology http://www.sciencedirect.com/science/article/pii/0921509389J 

142144, 1034 (2001), ISSN 0257-89 72 , 24 A. F. Wells, Structural Inorganic Chemistry, 5th ed. (Ox- 

jce:title^Proceedings of the 7th International Con- ford University Press, 1984) 

ference on Plasma Surface Engineeringj/ce:title^, 2j M. Kanoun and S. Goumri-Said, Physics Letters A 

http://www. sciencedirect. com/science/article/pii/S0257893BC2101091'B (2007), ISSN 0375-9601, 

N. Gordillo, R. Gonzalez-Arrabal, A. lvarez Herrero, and http://www.sciencedirect.com/science/article/pii/S037596010f 

F. Agull-Lpez, Journal of Physics D: Applied Physics 42 , 26 Z here, and in what follows, refers to the number of chem- 

165101 (2009), stacks. iop. org/ JPhysD/42/165101 ical formula units per unit cell. 

R. Cremer, M. Witthaut, D. Neuschtz, 27 W. B. Pearson, Handbook of Lattice Spacings & 

C. Trappe, M. Laurenzis, O. Winkler, and Structures of Metals & Alloys (International Series 
H. Kurz, Microchimica Acta 133, 299 (2000), of Monographs on Metals Physics & Physical Met- 
ISSN 0026-3672, 10.1007/s006040070109, allurgy) (Franklin Book Co, 1964) ISBN 0080090788, 

http : //dx . doi . org/ 10 . 1007/s006040070109 http: //www . amazon. com/Handbook- Structure s-International-Mon 

T. Maruyama and T. Morishita, Applied Physics Letters 28 R. Yu and X. F. Zhang, 

69, 890 (1996), http://dx.doi.Org/10.1063/l.117978 Phys. Rev. B 72, 054103 (Aug 2005), 

M. Asano, K. Umeda, and A. Tasaki, http://link.aps.org/doi/10.1103/PhysRevB.72.054103 

Japanese Journal of Applied Physics 29, 1985 (1990), 29 H. Shimizu, M. Shirai, and N. Suzuki, 

http://jjap.jsap.jp/link7JJAP/29/1985/ Journal of the Physical Society of Japan 66, 3147 

L. Maya (1993) pp. 203-208, (1997), http : //jpsj . ipap . jp/link?JPSJ/66/3147/ 

http : //www . scopus . com/ inward/record . url?eid=2- s2 . 0- 002^2,1.1 l&9%ar^rppclQ=40Mid^.02a7«fi45^0iaagb3Ea6a2B4bflb6f 7126 

D. yuan Wang, N. Nakamine, and Y. Hayashi, skowski, Angewandte Chemie International Edition 
Journal of Vacuum Science & Technology A: Vacuum, Surfaces, aiffl, Films 4365 (2006), ISSN 1521-3773, 
16, 2084 (1998), http : //link . aip . org/link/? JVA/16/2084/1 http : //dx . doi . org/10 . 1002/anie . 200600431 

C. Gallardo-Vega and W. de la 31 R. Yu, Q. Zhan, and X. F. Zhang, 

Cruz, Applied Surface Science 252, Applied Physics Letters 88, 051913 (2006), 

8001 (2006), ISSN 0169-4332, http://link.aip.Org/link/7APL/88/051913/l 

http : //www. sciencedirect . com/science/article/pii/S0169 3 <i3&E050148&&i Barth and L. Hedin, 

J. Pierson, Vacuum 66, 59 (2002), ISSN 0042-207X, Journal of Physics C: Solid State Physics 5, 1629 (Feb 

http : //www . sciencedirect . com/science/article/pii/SOO422O7802$),(MI2l5$: //iopscience . iop . org/0022-3719/5/13/012/ 

Z. Hou, Solid State Sciences 10, 33 M. Pant and A. Rajagopal, Solid State Communications 

1651 (2008), ISSN 1293-2558, 10, 1157 (1972), ISSN 0038-1098, 

http : //www. sciencedirect . com/science/art icle/pii/S129325&8egO006Sw . sciencedirect . com/ science/art icle/pii/0038 109872? 

S. Terada, H. Tanaka, and K. Kub- 34 G. Kr esse and J. Hafner, 

ota, Journal of Crystal Growth 94, Physical Review B 47, 558 (Jan 1993), 

567 (1989), ISSN 0022-0248, http://link.aps.org/doi/10.1103/PhysRevB.47.558 

http : //www. sciencedirect . com/science/article/pii/00220 3 l46S9900389 Kresse and J. Hafner, 

M. G. Moreno-Armenta, A. Martnez- Physical Review B 49, 14251 (May 1994), 

Ruiz, and N. Takeuchi, Solid State Sciences http://link.aps.org/doi/10.1103/PhysRevB.49.14251 



17 



G. Kresse and J. Furthm- ' M. Born and J. R. Oppenheimer, Annalen der Physik 84, 

Her, Computational Materials Science 457 (1927) 

6, 15 (1996), ISSN 0927-0256, 58 R. M. Martin, Electronic Structure, Baste Theory and 

http : //www. sciencedirect . com/science/art icle/pii/0927025$9&d&2QB(Methods (Cambridge University Press, 2004) c.f. 



G. Kresse and J. Furthmiiller, 

Physical Review B 54, 11169 (Oct 1996), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 54 . 11169 
J. Hafner, Journal of Computational Chemistry 
29, 2044 (2008), ISSN 1096-987X, 

http : //dx . doi . org/10 . 1002/ j cc . 21057 
G. Kresse and D. P. Joubert, 

Physical Review B 59, 1758 (Jan 1999), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 59 . 1758 
K. Capelle, Brazilian Journal of Physics 
36, 1318 (12 2006), ISSN 0103-9733, 



Appendix L: Numerical Methods 

R. P. Feynman, Physical Review 56, 340 (Aug 1939), 
http : //link . aps . org/doi/10 . 1103/PhysRev . 56 . 340 

G. Grimvall, Thermophysical Properties of Materials 
(North Holland, 1986) 

Eq. 5 results in a negative E co h- However, another con- 
vention with positive E co h is also common, where energy 
signs in Eq. 5 change. 

Y. Mishin, D. Farkas, M. J. Mehl, and D. A. Papa- 
constantopoulos, Physical Review B 59, 3393 (Feb 1999), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 59 . 3393 



http : //www . scielo . br/scielo . php?scr ipt=sci_arttext&pio 6 ^S&03^(S3gp06MD70M^taHm=is3md J 



Hohenberg 



Physical Review 



136, 



and 
B864 



W. 

(Nov 



Kohn, 
1964), 



http : //link . aps . org/doi/10 . 1103/PhysRev . 136 . B864 



W. Kohn 
Physical Review 



and L. J. Sham, 

140, A1133 (Nov 1965), 
http : //link . aps . org/doi/10 . 1103/PhysRev . 140 . A1133 
O. Gunnarsson and B. I. Lundqvist, 
Physical Review B 13, 4274 (May 1976), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 13 . 4274 
H. J. Monkhorst and J. D. Pack, 
Physical Review B 13, 5188 (Jun 1976), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 13 . 5188 
O. Jepson and O. Anderson, Solid State Communications 
9, 1763 (1971), ISSN 0038-1098, 

http : //www. sciencedirect . com/science/art icle/pii/OO381O980(H9iD,31EBmation 
G. Lehmann and M. Taut, physica status solidi (b) 
54, 469 (1972), ISSN 1521-3951, 

http : //dx . doi . org/10 . 1002/pssb . 2220540211 
P. E. Blochl, O. Jepsen, and O. K. Ander- 
sen, Physical Review B 49, 16223 (Jun 1994), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 49 . 16223 
M. Methfessel and A. T. Paxton, 
Physical Review B 40, 3616 (Aug 1989), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 40 . 3616 
A. D. Becke, Physical Review A 38, 3098 (Sep 1988), 
http : //link . aps . org/doi/10 . 1103/PhysRevA . 38 . 3098 
J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. 
Jackson, M. R. Pederson, D. J. Singh, and C. Fi- 
olhais, Physical Review B 46, 6671 (Sep 1992), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 46 . 6671 
J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. 
Jackson, M. R. Pederson, D. J. Singh, and C. Fi- 
olhais, Physical Review B 48, 4978 (Aug 1993), 
http : //link. aps . org/doi/10 . 1103/PhysRevB . 48 . 4978 . 2 
J. P. Perdew, K. Burke, and M. Ernzerhof, 
Physical Review Letters 77, 3865 (Oct 1996), 
http : //link . aps . org/doi/10 . 1103/PhysRevLett . 77 . 3865 
J. P. Perdew, K. Burke, and M. Ernzerhof, 
Physical Review Letters 78, 1396 (Feb 1997), 
http : //link . aps . org/doi/10 . 1103/PhysRevLett . 78 . 1396 
M. Ernzerhof and G. E. Scuseria, 
The Journal of Chemical Physics 110, 5029 (1999), 
http : //link . aip . org/link/? JCP/110/5029/1 
P. E. Blochl, Physical Review B 50, 17953 (Dec 1994), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 50 . 17953 
E. R. Davidson, Journal of Computational Physics 
17, 87 (1975), ISSN 0021-9991, 



Furthmuller, 

"Vasp the guide," (2011), available on-line at 
http: //cms .mpi .univie . ac . at/vasp/vasp/ . Last 
accessed October 2012. 

K. Lejaeghere, V. V. Speybroeck, G. V. Oost, 
and S. Cottenier, "Error bars for solid-state density- 
functional theory predictions: an overview by means 
of the ground-state elemental crystals," (August 2012), 
arXiv:1204.2733v2 [cond-mat.mtrl-sci] 
It is also well known that GGA may slightly lower the 
ground-state energy when a nonspherical ground-state 
density is allowed for (cf. Ref. [100] and Ref. 46 therein). 
N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Philadelphia: Saunders College, 1976) 
J. K. Dewhurst, "EOS version 1.2," (August 
of State (EOS) program for fit- 
ting energy- volume data (http://elk.sourceforge.net, 
http: / /exciting. sourceforge.net) . 

F. Birch, Physical Review 71, 809 (Jun 1947), 
http : //link . aps . org/doi/10 . 1103/PhysRev . 71 . 809 

D. A. Young, Phase Diagrams of the Elements (Uni- 
versity of California Press, 1991) ISBN 0520074831, 
http: //www . amazon. com/Phase-Diagrams-Element s-David-Young/d 
CRC Handbook of Chemistry and Physics, 82nd Edition 
(CRC Press, 2001) ISBN 0849304822 
J. Donohue, The Structures of the Elements 
(John Wiley & Sons Inc, 1974) ISBN 0471217883, 
http: //www . amazon. com/Structures-Elements- Jerry- Donohue/dp/' 
M. J. Mehl and D. A. Papaconstantopou- 
los, Phys. Rev. B 54, 4519 (Aug 1996), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 54 . 4519 
R. W. G. Wyckoff, The Structure of Crystals (The Chem- 
ical Catalog Co., New York, 1935) 

M. Sigalas, D. A. Papaconstantopoulos, and 
N. C. Bacalis, Phys. Rev. B 45, 5777 (Mar 1992), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 45 . 5777 
C. Kittel, Introduction to Solid State Physics, 5th edition, 
(John Wiley & Sons, New York, 1996) 
P. H. T. Philipsen and E. J. Baerends, 
Phys. Rev. B 54, 5326 (Aug 1996), 

http : //link . aps . org/doi/10 . 1103/PhysRevB . 54 . 5326 

E. Zarechnaya, N. Skorodumova, S. Simak, B. Johans- 
son, and E. Isaev, Computational Materials Science 
43, 522 (2008), ISSN 0927-0256, 
http: //www . sciencedirect . com/ science/art icle/pii/S092702560? 
S. Raju, E. Mohandas, and V. Raghunathan, J. Phys. 
Chem Solids 58, 1367 (1997) 

M. G. Moreno-Armenta, W. L. Prez, and N. Takeuchi, 



http : //www. sciencedirect. com/science/article/pii/0021999$SBa086SCb Sciences 9, 166 (2007), ISSN 1293-2558, 



18 



http://www.sciencedirect.com/science/article/pii/S129iiS5Sfi06fflJ2S58o/ Mineralogy, edited by J. W. An- 
J. Xiao, Y. Li, and A. Jiang, thony, R. A. Bideaux, K. W. Bladh, and M. C. 

Journal of Materials Science & Technology Nichols (Mineralogical Society of America, Chan- 

27, 403 (2011), ISSN 1005-0302, tilly, VA 20151-1110, USA) available on-line at 

http : //www. sciencedirect . com/science/art icle/pii/SlOO5O3B£tp600826 .handbookof mineralogy . org/ 



W. Yu, L. Li, and C. Jin, Journal of Materials Science 40, 
4661 (2005), ISSN 0022-2461, 10.1007/sl0853-005-0638- 
3, http : //dx . doi . org/10 . 1007/sl0853-005-0638-3 
J. G. Niu, W. Gao, X. P. Dong, L. Guan, and 
F. Xie, Advanced Materials Research (Volumes 150-151) 
Advances in Composites, 1290 (October 2010), 
www . scientific .net 

T. Maruyama and T. Morishita, 

Journal of Applied Physics 78, 4104 (1995), 
http : //dx . doi . org/10 . 1063/1 . 359868 
L. Maya, Journal of Vacuum Science & Technology A 
11, 604 (May 1993), http://dx.doi.Org/10.1116/l.578778 
N. Kanoun-Bouayed, M. BenaliKanoun, and S. Goumri- 
Said, Central European Journal of Physics 9, 205 
(2011), ISSN 1895-1082, 10.2478/sll534-010-0063-3, 
http : //dx . doi . org/10 . 2478/sll534-010-0063-3 
W. Hong-Bo and X. De-Sheng, 

Chinese Physics Letters 21, 1612 (2004), 

http : //iopscience . iop . org/0256-307X/21/8/056 
In Table I, our computed properties of the elemental Cu 
are compared with experiment and with previous calcu- 
lations as well. This may benchmark the accuracy of the 
rest of our calculations. 

The electronic structure of D0g , DO2 and RI1F3 phases are 
discussed in Sec. V, while the rest are not shown here. 
In fact, the accuracy of the approximate XC functional 
(PBE and others) does not really allow us to make a dis- 
tinction among these. 

M. Llunell, S. Alvarez, and P. Alemany, 
Journal of the Chemical Society, Dalton Transactions, 
1195(1998), http : //dx . doi . org/10 . 1039/A708200B 
M. F. C. Ladd, Introduction to Physical Chemistry 
(Cambridge University Press, 1998) ISBN 0521578817, 



M. Fuchs, J. L. F. Da Silva, C. Stampfi, J. Neugebauer, 
and M. Schemer, Phys. Rev. B 65, 245212 (Jun 2002), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 65 . 245212 

Table of interatomic distances and configuration in 
molecules and ions, edited by L. E. Sutton (The Chemical 
Society, London, UK, 1958) 

J. A. Montoya, A. D. Hernandez, C. San- 
loup, E. Gregoryanz, and S. Scandolo, 
Applied Physics Letters 90, 011909 (2007), 
http : //link . aip . org/link/?APL/90/011909/ 1 

A. F. Young, J. A. Montoya, C. Sanloup, 
M. Lazzeri, E. Gregoryanz, and S. Scan- 
dolo, Phys. Rev. B 73, 153102 (Apr 2006), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 73 . 153102 
C. Stampfi and C. G. Van de Walle, 
Phys. Rev. B 59, 5521 (Feb 1999), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 59 . 5521 
S. Narasimhan and S. de Gironcoli, 
Phys. Rev. B 65, 064302 (Jan 2002), 
http : //link . aps . org/doi/10 . 1 103/PhysRevB . 65 . 064302 

B. Grabowski, T. Hickel, and J. Neuge- 
bauer, Phys. Rev. B 76, 024309 (Jul 2007), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 76 . 024309 

C. J. Bradley and A. P. Cracknell, The Mathematical 
Theory of Symmetry in Solids: Representation Theory 
for Point Groups and Space Groups (Oxford: Clarendon 
Press, 1972) 

M. Gajdos and K. Hummer and G. Kresse 
and J. Furthmiiller and F. Bechstedt, 
Physical Review B 73, 045112 (Jan 2006), 
http : //link . aps . org/doi/10 . 1103/PhysRevB . 73 . 045112 
W. G. Aulbur, L. Jonsson, and J. W. 
Wilkins (Academic Press, 1999) pp. 1 218, 



http : //www. amazon. com/Introduct ion- Physical- Chemistry-MalitepLai6^i^dp/3ia2^i£^iBilryect . com/ science/art icle/pii/S008 119470? 



N. Wang, W.-Y. Yu, B.-Y. Tang, L.-M. 
Peng, and W.-J. Ding, Journal of Physics 
D: Applied Physics 41, 195408 (2008), 
http : //stacks . iop . org/0022-3727/41/i=19/a=195408 
C. Stampfi and A. J. Freeman, 
Physical Review B 67, 064108 (Feb 2003), 
http : //link. aps . org/doi/10 . 1103/PhysRevB . 67 . 064108 
P. Atkins and J. D. Paula, Atkins' Physical Chem- 
istry (Oxford University Press, 2009) ISBN 0199543372, 



J. Kohanoff, Electronic Structure Calculations for Solids 
and Molecules : Theory and Computational Methods 
(Cambridge University Press; Cambridge, 2006) 
J. Harl, The Linear Response Function in Density Func- 
tional Theory: Optical Spectra and Improved Description 
of the Electron Correlation, Ph.D. thesis, University of 
Vienna (2008), http://othes .univie . ac . at/2622/ 
The physical meaning of the RPA is that electrons are 
considered to respond to the total (external plus induced) 



http : //www. amazon. com/ Atkins-Physical- Chemistry- Julio- Pe"fariiViipy(Qi9£(£Aa3y5 

C. Colinet, Intermetallics 11, 1095 (2003), ISSN 113 Recall that the optical region (visible spectrum) is about 
0966-9795, special Issue Dedicated to the International (390 ~ 750) nm which corresponds to (3.183 ~ 1.655) eV. 

Meeting on Thermodynamics of Alloys - TOFA 2002, 114 G. Onida, L. Reining, and A. Ru- 
http : //www. sciencedirect. com/science/article/pii/S096697&fi£(300ia^«. Mod. Phys. 74, 601 (Jun 2002), 



H. W. Hugosson, P. Korzhavyi, U. Jans- 
son, B. Johansson, and O. Eriksson, 
Physical Review B 63, 165116 (Apr 2001), 
http : //link. aps . org/doi/10 . 1103/PhysRevB . 63 . 1651 16 
If Ecoh is used with a positive sign convention, i.e. negative 
of Eqs. 4 and 5, then signs in Eq. 11 must be reversed. 
M. D. Graef and M. E. McHenry, Structure of ma- 
terials : An Introduction to Crystallography, Diffrac- 
tion and Symmetry (Cambridge University Press, 2007) 
http : //www . Cambridge . org/9780521651516 



http : //link . aps . org/doi/10 . 1103/RevModPhys . 74 . 601 

A. D. Becke, The Journal of Chemical Physics 98, 5648 
(1993), http : //link . aip . org/link/? JCP/98/5648/1 
K. N. Kudin, G. E. Scuseria, and R. L. Mar- 
tin, Phys. Rev. Lett. 89, 266402 (Dec 2002), 
http : //link . aps . org/doi/10 . 1103/PhysRevLett . 89 . 266402 
E. N. Brothers, A. F. Izmaylov, J. O. Nor- 
mand, V. Barone, and G. E. Scuseria, 
The Journal of Chemical Physics 129, 011102 (2008), 
http : //link . aip . org/link/? JCP/129/011102/1 



19 



W.-G. Yin, D. Volja, and W. Ku, 
Phys. Rev. Lett. 96, 116405 (Mar 2006), 
http : //link . aps . org/doi/10 . 1103/PhysRevLett . 96 . 1 16405 

C. Franchini, R. Kovacik, M. Marsman, S. S. 
Murthy, J. He, C. Ederer, and G. Kresse, "Maxi- 
mally localized wannier functions in LaMn03 within 



PBE+U, hybrid functionals, and partially self-consistent 
GW: an efficient route to construct ab-initio tight- 
binding parameters for eg perovskites," (May 2012), 
arXiv: 1 1 1 1 . 1 528v2 [cond-mat . str-el] 



